Tuning the fragility of a glassforming liquid by curving space 
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We investigate the influence of space curvature, and of the associated "frustration" , on the dynam- 
ics of a model glassformer: a monatomic liquid on the hyperbolic plane. We find that the system's 
fragility, i.e. the sensitivity of the relaxation time to temperature changes, increases as one decreases 
the frustration. As a result, curving space provides a way to tune fragility and make it as large as 
wanted. We also show that the nature of the emerging "dynamic heterogeneities", another distinc- 
tive feature of slowly relaxing systems, is directly connected to the presence of frustration-induced 
topological defects. 



Among the many anomalous properties associated 
with glass formation, "fragility" is one that has attracted 
much attention [U El O SHU H [7] . Large fragility, i.e. 
large deviation of the temperature dependence of the 
viscosity and of the structural relaxation time from an 
Arrhenius behavior, is usually taken as the signature of 
a collective phenomenon that grows as temperature de- 
creases. This is certainly one incentive for the continuing 
search for a theory of the glass transition [2j [4j |5j |6j [7] . 
Yet, the absence of a simple glassforming liquid model in 
which one can control the degree of fragility, hence the ex- 
tent to which collective behavior develops, has hindered 
progress on developing and testing candidate theories. 

Since the early work of Frank [8 , a promising line of re- 
search on supercooled liquids and the glass transition has 
relied on the concept of "geometric frustration" [7J [9j [10] . 
Frustration in this context can be defined as an incom- 
patibility between extension of the local order preferred 
in a liquid and tiling of the whole space. The paradigm 
is the icosahedral order in metallic liquids and glasses, 
which although locally favored cannot tile space due to 
topological reasons [5]. Frustration of the icosahedral or- 
der, however, can be suppressed by leaving the Euclidean 
world and curving space [9j [10]. In a series of insightful 
articles[10j [TTJ [12], Nelson and collaborators have pro- 
posed a simpler two-dimensional (2D) analog: by placing 
a liquid of disks on a 2D manifold of constant negative 
curvature (the hyperbolic plane) , the local hexagonal or- 
der that can tile the ordinary Euclidean plane is now 
frustrated in a way which mimics by many aspects the 
frustration of icosahedral order in 3D Euclidean space. 
The model of a monatomic liquid on the hyperbolic plane 
therefore offers the opportunity to investigate, at a mi- 
croscopic level, the influence of the degree of frustration, 
here controlled by the curvature, on the slowing down of 
the relaxation associated with glass formation. 

We present the results of the first computer simula- 
tion of the dynamics of a liquid in curved hyperbolic 
space. The hyperbolic plane H 2 , also called pseudo- 
sphere or Bolyai-Lobatchevski plane, is a Riemannian 
surface of constant negative curvature [I3j [14] • Contrary 
to a sphere, which is a surface of constant positive cur- 



vature, H 2 is infinite: this allows one to envisage the 
thermodynamic limit at constant curvature. However, 
H 2 cannot be embedded as a whole in the 3D Euclidean 
space and "models", i.e. projections, must be used for 
its vizualization. The hyperbolic metric is often given in 
polar coordinates (r, 0), namely, 



ds 2 = dr 2 



sinh(ftr) 



d</> 2 
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which makes apparent the connection with the more fa- 
miliar metric of the sphere S 2 that is obtained by replac- 
ing the parameter n by in. The Gaussian curvature of 
H 2 is — n 2 < 0; n therefore measures the deviation from 
fiat space and n~ x can be taken as an intrinsic frustration 



length. We consider the Poincare disk model (Fig. 1(a)) 



that maps the whole infinite space H 2 onto the open disk 
of radius unity. This projection (r' = tanh(^r/2), <j> f = (j>) 
is conformal, i.e. it preserves the angles, but is not iso- 
metric: the Euclidean distance between two points of the 
disk separated by a given distance in H 2 shrinks to zero 
when the points approach the disk perimeter. 

To carry out a Molecular Dynamics (MD) simulation 
of particles on the hyperbolic plane, a number of serious 
methodological problems have to be resolved, which we 
only briefly allude to. Once the model is properly defined, 
the two main ingredients in any MD simulation are the 
algorithm to solve the Newton equations of motion and 
the boundary conditions, usually chosen as periodic in or- 
der to more rapidly converge to the thermodynamic limit 
corresponding to the experimental situation. Among the 
peculiarities one encounters when leaving flat space to 
consider curved manifolds such as H 2 is the absence of a 
global definition of parallel vectors. We handle this and 
generalize the standard MD algorithm to the hyperbolic 
plane by using a method detailed in a forthcoming publi- 
cation. Even more delicate is the question of the periodic 
boundary conditions (pbc's). Due to the hyperbolic na- 
ture of the metric, the contribution of the boundary of 
any finite system is always of the same order of magni- 
tude as that of the bulk of the system. Implementing 
proper pbc's is therefore crucial. Again, one must ac- 
count for the specificities of hyperbolic space: first, an 
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infinite number of regular tilings of H 2 are possible, and 
second, the area of the elementary cell of a given tiling 
is fixed by the curvature (see below) [T3l [14] . As a conse- 
quence of the latter property, studying finite-size effects 
at constant curvature requires to change the boundary 
condition. Building on our earlier work[15 , we have im- 
plemented two different pbc's: an octagonal pbc (Fig. 
1(a) ) and a pbc with a larger unit cell formed by a reg- 
ular 14-gon with a specific pairing of the edges. 

We have studied a monatomic liquid model in which 
atoms pairwise interact via the standard truncated 
Lennard- Jones potential v(r) = 4e ((a/r) 12 — (a/r) 6 ) 
(with a cut-off at 2.5a), where the distance r is defined 
with the hyperbolic metric. The control parameters are 
the temperature T, the density p, and the frustration as- 
sociated with space curvature and characterized by the 
dimensionless parameter kg. As already stressed, there 
are a number of geometrical and topological constraints 
associated with the pbc. For instance, the area of the 
unit cell is A = AitK~ 2 (g — 1), where g is the genus of 
the compact manifold associated with the pbc and is 
equal to 2 for the octagon and 3 for the 14-gon. For 
a given density, the number N of atoms is thus fixed by 
the curvature and the pbc: in the present work, N typi- 
cally varies between 300 and 30 000 as one decreases the 
frustration. We are interested by weak frustration for 
which the local order remains hexagonal as in the Eu- 
clidean plane. Indeed, for large enough frustration, the 
preferred arrangement of atoms around a central one is 
no longer a hexagon, but a heptagon, with unfrustrated 
extension of the hept agonal order to the whole space; on 
further increasing the frustration, one encounters locally 
preferred arrangements formed by polygons with an in- 
creasingly larger number of sides [12] . (Disk packings for 
large negative curvature have also been recently consid- 
ered in [16].) Simulations are performed for kg spanning 
one order of magnitude from 0.02 to 0.2. In addition 
to computing usual static quantities, e.g. the pair cor- 
relation function, we have carried out a direct analysis 
of the topological defects, which will be discussed below, 
and we have monitored several dynamic observables char- 
acterizing the motion of the atoms. From the distance 
travelled by any atom j between two times t' and t f + t, 
dj{t' , t' + t), we compute the hyperbolic generalization of 
the self intermediate scattering function, 
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F s (k,t) = -^(P_i + ^(cosh(^-(0,t)))), 



(2) 



where P_ 



is a Legendre function of first kind (such 



that Eq. (2) reduces to the conventional spatial Fourier 
transform in the Euclidean limit [17]). 

As a benchmark, we have first considered the Eu- 
clidean case (k = 0). We find that, irrespective of the 
cooling rate, the liquid orders in an hexagonal structure 
at a temperature T*(pa 2 ) (e.g. T*(pa 2 = 0.85) ~ 0.75 



in Lennard- Jones units). No glass formation is therefore 
possible, as anticipated. A very different behavior is ob- 
served when frustration is switched on by curving space. 
No ordering transition takes place (the transition at T* 
is thus "avoided" [7, ) and the liquid phase can now be 
cooled in equilibrium below T*. 

From F s (k,t) with k chosen near the maximum of the 
static structure factor (k a -1 ), we have extracted the 
translational relaxation time r, which is determined when 
F s (k,t) =0.1 (see Fig. |l(b)| ). An alternative definition of 
r is obtained from the fit of F s (k, t) (beyond the plateau) 
to a stretched exponential, e~^/ r ^. Up to a multiplica- 
tive constant, the two definitions of r give similar results 
for the T dependence. At low enough T, one reaches the 
limit of the computer resources and the liquid falls out of 
equilibrium to freeze in an amorphous solid, i.e. a glass. 

As shown in Fig. [2] one observes a striking pattern of 
variation of fragility with frustration. At high T above 
T*, the data show no dependence on curvature. This 
is easily understood by combining the fact that the re- 
laxation remains a local phenomenon and that locally 
H 2 appears flat for the atoms (recall that kg <C 1). 
A marked deviation from Arrhenius dependence, i.e. a 
super- Arrhenius behavior, is found below T*. The mag- 
nitude of this deviation unambiguously increases as the 
frustration parameter is reduced and one gets closer to 
the unfrustrated Euclidean case. As predicted by the 
frustration-limited domain theory [7], fragility therefore 
increases as frustration decreases. A rationale for this 
trend is that as the system gets closer to the avoided 
transition, the spatial correlations associated with frus- 
trated ordering grow larger; collective behavior thus oc- 
curs on longer length scales, which results in a more 
strongly super- Arrhenius dependence of the relaxation 
time and a larger fragility. Note that in line with this 
fragility pattern, the stretching exponent f3 is found to 
decrease with decreasing frustration: at T/T* c± 0.85, 
(3 = 0.64, 0.54, 0.50, 0.42 for kg = 0.2, 0.1, 0.05, 0.02. 

A crude heuristic argument suggests that the increase 
of fragility goes logarithmically with the inverse of the 
frustration. The idea is to compare the energy scales 
involved in the activation barriers for relaxation at low 
and high T and derive an estimate of the fragility through 
their ratio. At high T, the scale is provided by the in- 
teraction energy between atoms and is independent of 
curvature as seen from Fig. [2] On the other hand, one 
expects that the low-T dynamical behavior is controlled 
by the motion of the rare frustration-induced defects (see 
below). An estimate for the associated energy scale is 
obtained by considering the continuum approach valid 
at sufficiently large wavelength and low T. The energy 
of frustration-induced disclinations (in an otherwise hex- 
atic medium) is found to be a constant plus a term pro- 
portional to ln(l/ tanh(^(j/2))[T0], which for small frus- 
tration behaves as ln(l/(ft<r)). The ratio of the energy 
scales at low and high T therefore goes as the logarithm of 
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FIG. 1: (a) The projection of H 2 associated with the Poincare disk model: atoms are represented by disks whose diameter 
contracts as they move away from the center of the Poincare disk. The geodesies are arcs of circles perpendicular to the disk 
boundary. We also display the octagonal tiling (in which 8 octagons meet at each vertex) used in the pbc. (b) F s (k,t) versus 
t for kg — 0.2 and for different temperatures T/T* ranging from 3.9 to 0.1 (left to right): for T ~ T* and below, stretched 
exponential behavior, exp(— (t/r)^), is observed with (3 decreasing from 1 at the highest T to 0.54 at the lowest equilibrated 
T; for T < 0.35 T* the system is out of equilibrium and ages (last three curves). 



l/(ft<r), which is compatible with the variation of fragility 
obtained from the simulation data. This indicates that 
fragility can be made as large as wanted by taking the 
limit of vanishingly small curvature. 

Another canonical feature of slowly relaxing systems, 
glassforming liquids in particular, is the "heterogeneous" 
nature of the dynamics [H El EDI EQ. This Phe- 
nomenon is easily detected by following the particle tra- 
jectories for a given period of time, as shown in Figs. 
3 (a),(b). Topological defects [22] and/or medium-range 
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FIG. 2: Arrhenius plot of the translational relaxation time r 
versus T*/T for several curvatures, where T* is the (approx- 
imate) location of the ordering transition in the Euclidean 
plane. The particle density is pa 2 = 0.85. Note the devi- 
ation from the simple Arrhenius dependence shown by the 
dashed line. As the frustration kg decreases, the deviation 
is stronger, i.e. the fragility gets larger. The thick blue line 
above T* corresponds to the Euclidean case. The continuous 
lines are fits to the VFT formula, r = r exp(L>T /(T - T )): 
as frustration decreases, D is found to decrease from 10.8 to 
4.6 while To/T* increases from 0.12 to 0.33, which again in- 
dicates an increase in fragility. 



ordering [23 have been suggested as playing a role in dis- 
ordered 2D phases exhibiting dynamic heterogeneities, 
and the present 2D monodisperse glassforming liquid of- 
fers a unique opportunity to investigate this point. 

The topological defects can be defined at a microscopic 
level by analyzing the local environment of each atom. To 
do so, we use a curved-space generalization of the Voronoi 
tessellation that provides an unambiguous means to as- 
sign a number of nearest neighbors to each atom. Most 
atoms have 6 neighbors, which corresponds to the hexag- 
onal local order. Defects, more specifically point "discli- 
nations" , are then associated with atoms with a coordina- 
tion number different from 6. Negative curvature forces 
in an irreducible number of disclinations of negative topo- 
logical charge (more than 6 neighbors), this number be- 
ing fixed by topological relations. In Figs. 3 (c),(d), we 
display the same atomic configurations as those plotted 
in Figs. 3 (a),(b) with a color code indicating the co- 
ordination number of each atom. At high T, there is a 
large density of defects and it is hardly possible to no- 
tice the imbalance in favor of negative disclinations. At 
low T on the other hand, the number of defects is small 
and one can clearly see the emergence of large domains 
of local 6-fold order coexisting with small localized de- 
fective regions. A closer inspection reveals that there are 
exactly 12 such regions consisting of a 7- fold disclination 
and attached short strings of little dipoles of 5-fold and 
7-fold disclinations forming "dislocations". Such strings 
have been dubbed "grain boundary scars" in the context 
of crystals on spherical surfaces [24] . (Note that the ir- 
reducible number of 7-fold disclinations is fixed by the 
pbc, here 12 (g — 1) with g = 2, but that the irreducible 
density of disclinations decreases, and the typical size of 
the locally ordered domains increases, as curvature de- 
creases.) Comparison of the two sets of figures clearly 
shows that the emerging heterogeneous character of the 
dynamics is directly linked to the topological defects: at 
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FIG. 3: (a),(b) Atomic trajectories followed for a time interval 
during which the average distance travelled by the atoms is 
roughly a at (a) T/T* = 2.4 and (b) T/T* = 0.52. Whereas 
at high T all atoms seem to move by a comparable amount 
and the dynamics is thus spatially homogeneous, a strikingly 
different picture is obtained at low T: most atoms hardly 
move or just rattle in the cage formed by their neighbors and 
mobility is concentrated in rare localized regions, illustrat- 
ing the spatial heterogeneity of the dynamics. (c),(d) Same 
atomic configurations as in (a) and (b) (respectively), with 
the color indicating the coordination number for each atom. 
Black: 6 neighbors (local hexagonal order), red: 7 neighbors 
(disclination of topological charge —tt/3 ), blue: 5 neighbors 
(disclination of charge +7r/3); at T/T* =2.4 (c), there are 
also defects with larger charges (orange: 8 neighbors, cyan: 
4 neighbors). Note the correspondence at low T between the 
rare localized defective regions in (d) and the regions of high 
mobility in (b). The frustration is kg — 0.05 and the particle 
density pa 2 = 0.85. 



low T, the regions of high mobility coincide with the 
vicinity of the intrinsic frustration-induced defects and 
their attached strings of dislocations whereas the regions 
of low mobility coincide with the hexagonal patches. We 
stress that the system is in a liquid phase even at low 
T and that all defects and atoms move over long enough 
time spans. Preliminary results on the 4-point space- 
time correlation function x±(fy obtained as the variance 
of the local relaxation associated with Fg(k. t)[25] indi- 
cate that the spatial correlations in the dynamics are 
maximum around a timescale of the order of the relax- 
ation time r and increase continuously as T decreases. 
The phenomenology is thus similar to that found in other 
glassformersOSl [RJ EDI HQ • 

The above results suggest the passage as T decreases 



from a local atomic dynamics to a collective relaxation 
controlled by the motion of topological defects, with an 
intermediate region that is determined by the proxim- 
ity to the avoided transition at T*. This intermediate 
region becomes more important as frustration decreases 
and one expects that in this regime, growth of static spa- 
tial correlations, super- Arrhenius behavior and extension 
of the dynamic heterogeneities all go together. How- 
ever, the extent of static spatial correlation saturates to a 
value given by the intrinsic frustration length n~ x (which 
fixes the average distance between the remaining intrin- 
sic disclinations) ; at low T and over distances beyond 
ft -1 , the slowing down of relaxation should therefore be 
controlled by the rare intrinsic 7-fold disclinations. One 
may speculate that in this regime the growing dynamic 
correlations (as extracted from x±(P)) decouple from the 
static spatial correlations and reflect longer-range cor- 
relation among the mobilty regions associated with the 
residual defects. Work is now in progress to investigate 
this potential low-T decoupling phenomenon. 
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